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Abstract 

We study heteroclinic networks in R^, made of a certain type of 
simple robust heteroclinic cycle. In simple cycles all the connections 
are of saddle-sink type in two-dimensional fixed-point spaces. We 
show that there exist only very few ways to join such cycles together 
in a network and provide the list of all possible such networks in R'^. 

The networks involving simple heteroclinic cycles of type A are new 
in the literatnre and we describe the stability of the cycles in these 
networks: while the geometry of type A and type B networks is very 
similar, stability distingnishes them clearly. 
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1 Introduction 

The study of robust heteroclinic cycles and networks is well-established as an 
interesting subject in the scientific community. It was approached initially 
from the point of view of the existence of such heteroclinic object^ as in 
Guckenheimer and Holmes and dos Reis [28]. Then, the study of the 
asymptotic stability of cycles led to results in Krupa and Melbourne [191 
[2U] . When addressing the stability of cycles in networks, it is clear that 

^Long-term address: Fachbereich Mathematik, Universitat Hamburg, BundesstraBe 55, 
20146 Hamburg, Germany 

^We point out that the purpose of the references in this section is to illustrate rather 
than to exhaust the existing bibliography on the subject. 
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no individual cycle can be asymptotically stable and intermediate notions 
of stability appeared in Melbourne [22], Brannath |1], Kirk and Silber |18j . 
Driesse and Homburg [9] and Podvigina and Asliwin [25] . Of these essential 
asymptotic stability (e.a.s.jl is the strongest and the one we use. 

In another spirit, an interest was taken in the dynamics near networks 
of various degrees of complexity, see e.g. Homburg and Knobloch [T5] or 
Postlethwaite and Dawes [2Z]. Clearly, this must be done on a case-by-case 
basis although an efficient tool seems to be the use of symmetry to study 
a quotient network instead of the original one. Examples may be found in 
Aguiar et al. [2], Aguiar and Castro [T], and Castro et al. [5]. 

The manifest interest in, and complexity of, the dynamics observed near 
networks provides a good reason for a systematic treatment of such objects, 
particularly in low dimensions where the analysis is tractable. At this point, 
we must distinguish between homoclinic and heteroclinic cycles: the former 
exhibit connections of a node to itself (or to another node in the same group 
orbit). Homoclinic cycles were systematically addressed by Sottocornola [31] 
and Podvigina [21] (see also Homburg et al. m), while networks involving 
homoclinic cycles appear in the results of Driesse and Homburg n and Pod¬ 
vigina and Chossat [26] . We focus on a particular type of heteroclinic cycle, 
in the context of symmetry, made of what we call elementary building blocks, 
see section [21 These cycles are such that the equilibria lie on the coordinate 
axes which are hxed-point spaces, the connections occur in coordinate planes 
(also hxed-point spaces) and all the eigenvalues are real. See Figures [H and 
12] for an illustration. 

A systematic study of heteroclinic networks in low dimension is useful 
in that it provides simple examples of interesting dynamics and a gallery 
of case studies. The study of bifurcations from heteroclinic networks, for 
instance, can greatly beneht from such a gallery. See Kirk et al. [I?] for 
an interesting example of such a phenomenon. A hrst step towards such a 
systematic approach may be found in Castro and Lohse [6] . The construction 
of heteroclinic networks in the context of coupled cell systems has recently 
been addressed by Ashwin and Postlethwaite [3j and by Field m, with 
a different purpose. Given a heteroclinic network these authors hud ways 
in which a coupled cell system may be constructed so that its dynamics 
exhibit the prescribed network. The present article provides a list of certain 
simple networks in made of simple cycles in the sense of 123. satisfying 
Assumption A at the beginning of section [3l Whether and how these can be 
realized in the context of coupled cell systems is the object of [3] and P. 
The simplex and cylinder methods of [31 Section 2] are of particular interest 

^In |25] this is called predominant asymptotic stability. 
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since they ensure the networks have most of the desired properties. 

This article’s contribution consists of a complete list of a special type of 
heteroclinic network in dimension 4, together with the study of the stability 
properties of each individual cycle. In particular, heteroclinic networks con¬ 
sisting of cycles of type A, so far absent from the literature, are considered. 
Although the geometry of networks of type A is similar to that of networks 
of type B, the stability properties of cycles are very different. 

Stability of individual cycles for type A networks is very constrained. In 
fact, for the particular networks we consider, all cycles except at most one 
attract almost no points in their neighbourhood. For an e.a.s. network of 
type A consisting of two cycles, while one of the cycles attracts virtually 
nothing in its neighbourhood, the other cycle attracts initial conditions near 
either cycle in the network. When the e.a.s. network is made of three cycles, 
two of these attract almost nothing in their respective neighbourhoods while 
the third cycle attracts initial conditions near any connection in the network. 
Which cycle has the attracting properties depends on the linearization at the 
nodes of the network. This is reminiscent of the type B network in the set-up 
studied by Kirk and Silber [IS] but very different from many cases addressed 
in [6]. 

From the point of view of applications the interest in heteroclinic cycles 
and networks manifests itself in subjects such as neuroscience, geophysics, 
game theory and populations dynamics. Illustrations in the literature appear 
in the work of Chossat and Krupa [7|, Rodrigues [22], Aguiar and Castro [T] 
and Hofbauer and Sigmund [T3|, respectively. 

By restricting to simple networks satisfying Assumption A in we ex¬ 
clude networks of potential interest but beyond the scope of this paper. These 
include those where nodes have complex eigenvalues, see for instance Ro¬ 
drigues and Labouriau [30] • Allowing for non-simple (again with respect to 
the dehnition in [20|) cycles there are networks with more than one equilib¬ 
rium on the same axis and on the same side of the origin. These can be ob¬ 
tained through the cylinder method of [3|. Another type of non-simple cycle 
occurs when the invariant space containing a connection is not a hxed-point 
space. Such examples appear for instance in replicator dynamics, see [TB] . 
In this case the existence of invariant subspaces is intrinsic to the dynamics 
rather than a consequence of symmetry. We can also obtain non-simple cy¬ 
cles by allowing equilibria outside the coordinate axes. Dynamics near such 
a network can be found in Kirk et ah |16j . 

The next section provides the necessary definitions and results for the 
remaining ones. Section [B] establishes simple results which lead to a complete 
list of a certain type of heteroclinic networks in R^. Examples of some of 
these are known in the literature. We provide the construction of the new 
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possibilities in an appendix. In Section 01 we resort to the calculations of 
Podvigina and Ashwin [2S] to establish the stability indices for each cycle in 
the network. Joining two or more cycles in a network places restrictions on 
the possible values given in |25] . Section [5] concludes. 

2 Preliminary results 

Our concern is with vector fields in described by a set of differential 
equations y = f{y), where / is P-equivariant for some finite group P C 0(4), 
that is, 

fi'y-y) = 7 -/( 2 /), V 7 G P V 1 / e 

A heteroclinic cycle consists of equilibria, also called nodes, i = 1,... ,m, 
together with trajectories which connect them: 

[e* ^ 6+i] c n iy*(6+i) ^0 (* = 1,...,m; = 6). 

If G r .^1 for all z = 2,... ,m, then we say the cycle is homoclinic. When 
referring to a cycle as heteroclinic we always assume it is not homoclinic. 
Such a cycle connects at least two nodes in different group orbits. It is clear 
that group orbits of equilibria in a cycle always appear in the group orbit of 
the whole cycle. See [23l subsection 2.5] for an interesting discussion on the 
definition of a heteroclinic cycle. 

Using the notion of a building block from [26] we describe the cycles we 
study as follows. A building block, see [2S1 P- 901], is a sequence of equilibria 
and connections [^i —?■ ... —)■ ^m+i] where ^m+i = 7^i for some 7 G P and 
such that ^ r .^1 for alH = 2,..., m. We say a building block is elementary 
if no two of the nodes ^ 1 ,..., ^^ belong to the same group orbit and 7 = Id. 

A hnite connected union of heteroclinic cycles is called a heteroclinic 
network. 

We assume that each connection [^i —)■ .^j+i] is of saddle-sink type in an 
invariant subspace in order to ensure robustness of the cycle. 

Heteroclinic cycles are classified as simple if the connections between con¬ 
secutive equilibria are contained in a two-dimensional subspace. We use the 
definition of Krupa and Melbourne [201 P- 1181]: let C P be an isotropy 
subgroup and let Pj = Fix(Ej). Assume that for all j = 1 ,... ,m the con¬ 
nection ^j+i] is a saddle-sink connection in Pj. Write Lj = Pj_i fl Pj. 

A robust heteroclinic cycle X C R^\{0} is simple if 

(i) dimPj = 2 for each j; 

(ii) X intersects each connected component of Lj\{0} in at most one point. 
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In most of the literature it seems to have been silently assumed that 
for simple cycles the linearization df{^j) has no double eigenvalues. We 
adopt this further assumption in referring to simple cycles and define simple 
networks accordingly. 

Definition 2.1. A heteroclinic network is called simple if 

(i) it is made of simple cycles, 

(ii) it intersects each connected component of Lj\{0} in at most one point. 

Simple cycles have been classihed into types A, B, and C, both in the 
context of bifurcation of cycles (see Chossat et al. 0) and in the context 
of their stability (see Krupa and Melbourne 121 and Podvigina and Ashwin 
pS]b A further type, Z, appears in Podvigina [23] to include some cycles of 
types B and C. We use the original classihcation into types A, B and C, 
which we reproduce here from 121. 

Definition 2.2 (Definition 3.2 in Krupa and Melbourne [2D|)- Let X C 
be a simple robust heteroclinic cycle. 

(i) X is of type A if Ej = Z 2 for all j. 

(ii) X is of type B if there is a hxed-point subspace Q with dimQ = 3, such 
that X G Q. 

(iii) X is of type C if it is neither of type A nor of type B. 

We use subscripts to indicate the number of distinct group orbits of equi¬ 
libria in a cycle, and superscripts to denote whether —Id is an element of P 
(—) or not (-I-). For example, a B^ cycle has three (group orbits of) equilibria 
and —Id G P. It follows from the dehnition (see also [201 Corollary 3.5]) that 
the vector helds supporting cycles of type A are equivariant under symmetry 
groups that do not possess a reflection, whereas those for cycles of types B 
or C do possess a reflection. In this context, reflection means reflection in a 
hyperplane. 

The asymptotic stability of the cycles, dependent on the eigenvalues of 
the vector field at each equilibrium, has been studied by Krupa and Mel¬ 
bourne na^. Joining cycles in a network prevents each cycle from being 
asymptotically stable, calling for intermediate notions of stability. These 
have been introduced by Melbourne Brannath |1], Kirk and Silber |18] . 
More recently, Podvigina and Ashwin [25] and Podvigina [23] have revisited 
these notions. 
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As in Podvigina and Ashwin [23], we denote by B^{X) an ^-neighbourhood 
of a (compact, invariant) set X C M"’. We write B{X) for the basin of at¬ 
traction of X, i.e. the set of points a: G M” with uj{x) C X. For 5 > 0 the 
5-local basin of attraction is := {x ^ B{X) \ (ptix) G Bs{X)yt > 0}, 

where 0t(.) is the flow generated by the system of equations. 

The following is the strongest intermediate notion of stability. We denote 
Lebesgue measure by £(.). 

Definition 2.3 (Dehnition 1.2 in Brannath |1]). A compact invariant set X 
is called essentially asymptotically stable (e.a.s.) if it is asymptotically stable 
relative to a set iV C M" with the property that 


t(BJX)) 


1 . 


In [25] the concept of stability index is introduced. It provides a means of 
quantifying the attractiveness of a compact, invariant set X. See Dehnition 
5 and section 2.3 in Podvigina and Ashwin Ba- 


Definition 2.4. For x ^ X and £, 5 > 0 dehne 


i{B,{x)nB{X)) ^ i{B,{x)nBs{X)) 

imx)) ’ • mix)) 


Then the stability index at x with respect to X is dehned to be 


aix) := cT+(x) — cT-(x), 


where 


a (x) := lim 

nn(S,(x))l 

, c, (x) := lim 

pn(l-S,(x))l 


[ ln(£) J 

£— >0 

[ ln(£) J 


The convention that (J-ix) = oo if = 0 for some e > 0 and o'+(x) = oo 

if Se(x) = 1 is introduced. Therefore, cr(x) G [—oo, oo]. In the same way the 
local stability index at x G X is dehned to be 


^loc(^) • *^loc,-|-(^) *^loc,—(^)) 


with 


o'ioc-(a;) 


lim lim 

S —^0 €—^0 


ln(Sg,^(x)) 

ln(£) 


*^loc,+ ix) 


lim lim 

S —>^0 s —^0 


ln(l - S£,^(x)) 
ln(£) 


6 

















The stability index o'(a;) quantifies the local extent (at x G X) of the 
basin of attraction of X. If o'(x) > 0, then in a small neighbourhood of x an 
increasingly large portion of points is attracted to X. If on the other hand 
ct ( x ) < 0, then the portion of such points goes to zero as the neighbourhood 
shrinks. 

Theorem 2.2 in |2S] establishes that both a{x) and crioc(x) are constant 
along trajectories. In order to characterize the attraction properties of a 
heteroclinic cycle in terms of the stability index we have to calculate only a 
finite number of indices. Podvigina and Ashwin 1251 denote the index along 
the trajectory leading to an equilibrium by cxj. In heteroclinic networks 
there may be more than one such connection, which is why in such cases 
we differ from this notation by writing cxjj for the index along the trajectory 
from to ^j. Moreover, Theorem 2.4 in [22] shows that the calculation of 
the indices can be simplified by restricting to a transverse section. 

Local stability indices are related to essential asymptotic stability in the 
following way. 

Theorem 2.5 (Theorem 3.1 [2T|). Let X C ML be a heteroclinic cycle or 
network with finitely many equilibria and connecting trajectories. Suppose 
that the local stability index ctioc(x) exists and is not equal to zero for all 
X G X. Then X is essentially asymptotically stable if and only if (7\oc{x) > 0 
along all connecting trajectories. 

Stability indices are always considered with respect to a set X. In the 
context of heteroclinic networks, this may be a single cycle or the entire 
network. We distinguish between these indices by referring to them as c- 
and n-indices, respectively. 

3 Construction of simple networks 

From now on we are concerned with simple networks in satisfying the 
following 

Assumption A: 

(i) All the cycles in the heteroclinic network have at least one elementary 
building block. 

(ii) All the equilibria in the heteroclinic network are on a coordinate axis. 

Assumption A(i) means that symmetric images of connections and equi¬ 
libria occur only as images of an entire cycle (or network). Assumption A(ii) 
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is naturally satisfied when the simplex and cylinder methods are used to 
generate simple cycles. Cycles of types B and C follow directly from these 
methods, cycles of type A can be obtained by adding some symmetry break¬ 
ing terms. 

In [HI Proposition 3.1], the authors have established some results con¬ 
cerning simple heteroclinic networks in involving cycles of types B and 
C. The results concern networks with at least one common connecting tra¬ 
jectory. We extend this study to address networks of type A. We also prove 
that having a common connecting trajectory is compulsory for networks sat¬ 
isfying Assumption A. The present section provides all the necessary results 
leading to this complete list. Note that Lemmas 13.11 to 13.31 do not require 
Assumption A. 

The next lemma completes Proposition 3.1 in [B] which deals with net¬ 
works made of only two cycles of types B and C. For the sake of completion, 
we include the results of [6] in the statement of Lemma 13.11 which therefore 
informs on all three types of network. 

Lemma 3.1. Cycles of type A can only be part of a network involving other 
cycles of the same type. Cycles of types B and C can he part of networks 
involving other cycles of the same type and egual number of eguilibria or by 
joining cycles of types Bf and Cf. 

Proof. According to the dehnition of the type of the cycles, the symmetry 
group of vector fields supporting cycles of type A does not possess any reflec¬ 
tions whereas that of vector fields supporting cycles of types B and C does. 
Therefore, cycles of type A are compatible only with cycles of the same type. 

Concerning cycles of types B and C, Krupa and Melbourne m section 
3.2] have provided a description of the symmetry groups of supporting vector 
helds. These are as follows 


symmetry group 

zi 

z^ 

Z 2 IK ^2 

z^ 

cycle of type 

^2^ 

^ 3 - 




concluding the proof. □ 

We now derive a few essential results that further limit the possibilities 
for heteroclinic networks in These are then put together to obtain a 
complete list of networks in Theorem 13.81 

Lemma 3.2. Consider a simple cycle with two nodes. These nodes belong 
to the same one-dimensional vector space. 

Proof. Let and ^2 be the nodes of a simple cycle [.^i ^2 ^ ■Ci]- Then 

[.^1 —)■ .^ 2 ] C P 2 and [.^2 ■Cl] P -Pi) where Pi, i = 1,2 is a two-dimensional 















fixed-point space. Also, G Pi fl P 2 = which is a one-dimensional 

vector space. □ 

Lemma 3.3. At each node of a simple network in there exist at most 
three connecting trajectories (up to symmetry). 

Proof. At each node of such a network there are four eigenvalues. Given 
that one of them is radial, there are three possible connecting directions 
corresponding to the remaining three eigenvalues. □ 

Lemma 3.4. Up to symmetry, a simple network in R^ satisfying Assumption 
A does not involve more than four eguilihria. 

Proof. Each connected component of Lj\{0} contains at most one equilib¬ 
rium in the network because the network is simple. So suppose for some i 
there are two equilibria G Li that are unrelated by symmetry. Then we 
have the following cases: 

(a) There is a connection [,^j —)■ within a plane Pij. Then Lj is not 

an invariant subspace, so it cannot contain equilibria and the overall 
number of nodes does not exceed four. 

(b) There are connections —)■ and within P^j. Then there 

exists 7 G T such that = 'yC,i, because at f,j the stable, unstable and 
transverse spaces are orthogonal complements. 

(c) There are connections [.^j —)■ and within P^j. Then and 

are in different connected components of Lj\{0}. It follows from [201 
Proposition 3.1] that, if P^j = Fix(E), then Z 2 C A^(S)/S and G T^,, 
contradicting our assumption. 

(d) There are connections [^i —?• and within planes P^ and 

Pik- Then the respective other connections (i.e. those leading to and 

fall into one of the cases above. 

Note that for none of the cases the directions of the connections matter. □ 

Corollary 3.5. In R^, the maximum number of connecting trajectories in a 
simple network satisfying assumption A is six (up to symmetry). 

Lemma 3.6. Two cycles in a simple network in R'^ satisfying Assumption 
A have at least one node in common. 
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Proof. Assume there are two cycles without a common node. By Lemma [23] 
each cycle has two equilibria, say —)■ ^2 [^f —)■ ^2 ■Cf]- 

Then there must be a connection between the two cycles, say ^f]. 

Because of Lemma 1331 for i = 1, 2 the are on one axis, say Li, and the 
are on another axis, say L 4 . The existence of connections between and ^2 
prevents the other axes from being hxed-point spaces. Therefore, the group 
orbit of Li is either itself or L 4 . 

If it is Li, then the group image of L 4 is also L 4 . Hence, T hxes Li and L 4 , 
so that neither is a hxed-point space, contradicting the dehnition of simple 
cycles. If it is L 4 , then T,^/^ — 1 there is only one cycle. □ 

Corollary 3.7. Two cycles in a simple network in satisfying Assumption 
A have at least one connecting trajectory in common. 

Proof. By Lemma [3.61 there is a common node fi. Suppose there is no com¬ 
mon connection to or from Then there must be two incoming and two 
outgoing trajectories from This contradicts Lemma [3.31 □ 

Theorem 3.8. In the following is the complete list of simple heteroclinic 
networks satisfying assumption A: 

• (^2,^2); (A3, A3); (A3,A4); (A3,A3,A4); 

Proof. Proposition 3.1 in [ 6 ] shows that the only networks made up of two 
cycles of type B 01 C are the ones listed here. Recall that the network 
is given in [B], as is the network (Rg , 6*4 ), while (R 3 , Rg ) is studied 
by Kirk and Silber (TB]. Note that the assumption of a common connecting 
trajectory in [ 6 l Proposition 3.1] is redundant given Corollary 13.71 Although 
the assumption on non-existence of critical elements outside the network in 
[HI Proposition 3.1] is stronger than Assumption A, we remark that the latter 
is sufficient for its proof. 

In order to show that we do not have any other networks involving types 
R and C, it follows from Lemma [2^ that constructing a network from a cycle 
requires the addition of a connection to an already existing node. Then, 

(i) adding a cycle to the (R^, Bf) network requires an additional connect¬ 
ing trajectory at one of the nodes, contradicting Lemma [3.31 

(ii) adding a cycle to the (Rj", Bf) or (Rj", C^) network, both of which al¬ 
ready have hve connecting trajectories, without contradicting Corollary 
13.51 leads to the (Rj", Rj", Cf) network. 
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The existence of the ,3^,0^^) network is established by Brannath [1]. 

From the previons resnlts we know that, other than the networks listed, 
only networks joining cycles of type A 2 to cycles of type either A 3 or A^ are 
possible. We begin by showing that there is an obstruction to the existence 
of such networks. 

Let Pij = {(xi,X 2 ,X 3 , 0 : 4 ) G I Xfc = 0 for k ^ {bj}}- By Lemma 
we may assume that the A 2 cycle is [.^1 —)■ ^2 —t .^ 1 ] where ^ 1,^2 £ Li. Also, 
we may assume that [^1 —)■ ^ 2 ] C P 12 and [^2 —^ ^ 1 ] C P 14 . From Corollary 
[33 assume that [^1 —)• ^ 2 ] C Pu is the common connecting trajectory. 

Then the A 3 cycle connects [^1 —)■ ^2 —^ ^3 —^ ^ 1 ] where we can choose 
[^2 ^ 3 ] to be in P 13 . Because and ^2 are on the same line, it follows 

that [^3 —)■ ^ 1 ] is also in P 13 . But then the connections in P 13 are not of 
saddle-sink type and therefore, are not robust. 

Concerning a connection to a cycle of type A 4 , note that since .^1 and ^2 
are on opposite sides of the origin, the connections between them in P 12 and 
Pi 4 must cross the X 2 - and a; 4 -axis, respectively. There can be no equilibria on 
these axes as otherwise they would be invariant, preventing the existence of 
the connections in P 12 and P 14 . Then the A 4 cycle must have two additional 
equilibria ^ 3,^4 on the a: 3 -axis and [^2 —>■ ^ 3 ], [^4 —^ ^ 1 ] C P 13 . Note that, 
for these cycles, the group orbit of the equilibria is non-trivial and 7 .^* must 
belong to either of the axes already containing equilibria. This contradicts 
(ii) in the definition of simple networks. 

In order to show that the listed networks exist, we provide symmetry 
groups with the required fixed-point spaces for the existence of the networks 
of type A. In appendix A we construct vector fields supporting the net¬ 
works of type A. A systematic approach to the construction of simple cycles, 
including homoclinic cycles, may be found in [ 20 ] where the quaternionic 
presentation is used as an alternative description of the symmetry groups. 

Since the symmetry group supporting a heteroclinic cycle of type A has 
no reflections but still has subgroups isomorphic to Z 2 , we construct a group 
generated by rotations by vr on coordinate planes. This is multiplication by 
— 1 of two of the four coordinates of and is isomorphic to Z 2 = (—Id). 
We look for groups generated by elements Kij such that Pij =Fix((Kjj)). For 
instance, we have 

K 12 = (Id, 

where is a rotation by tt on the plane 

^34 = {(0,0,0:3,0:4)}. 

Then ^ 12 ( 0 : 1 , 0 : 2 , 0 : 3 , 0 : 4 ) = {xi,X 2 ,-X 3 ,-X 4 ). 
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We construct the groups so as to guarantee the existence of the invariant 
planes for the heteroclinic connections. These are 

• the (^12,^2) network has one cycle with connections [^1 —)■ ^2 ^1] in 

P 12 U Pi3 and the other with connections [^1 —)■ .^2 —t ■Ci] -P12 U P14; 

• the (A 3 , A 3 ) network has one cycle with connections [^1 —t ^2 —t ^3 ^ 

^ 1 ] in P 12 UP 23 UP 13 and the other with connections [^1 —t ^2 —t ^4 ^ 1 ] 

in P 12 U P 24 U P 14 ; 

• the (A3, A4) network has the A3 cycle with connections [^1 ^ ^2 ^ 

^3 ^ 1 ] in P 12 U P 23 U Pi 3 and the A4 cycle with connections [^1 

^2 ^3 —)■ ^4 ^ ^1] in P12 U P23 U P34 U P14; 

• the (A3,A3,A4) has, additionally to the previous, the A3 cycle with 
connections [^1 —t ^2 —t ^4 —)• ^1] in P12 U P24 U P14. 

The groups for the respective networks are as follows: 

• the (A2, A2) network exists under the action of T = (ki2, ^13); 

• the (A3, A3)-, (A3, A4)- and (A3, A3, A4) networks exist under the action 
of T = (ki2, ksa)- 

We have ku = K12 o ^43, so T = (^12, K13) has P12, P13 and P14 as invariant 
planes. Since K23 = K13 o ^34 and ^24 = o %4, the group T = (^12, K13, K34) 
has as invariant planes all the coordinate planes. □ 

Considering the symmetry groups in the construction in the appendix, we 
hnd that —Id = ^12^34 G T for the type A cycles with three and four nodes. 
For A cycles with two nodes —Id cannot be an element of T since and .^2 
are on the same axis and —Id maps to —^j. In view of this, we write A^, 
Aj and Aj from now on0 

The networks made of cycles with the same number of nodes of type A or 
B only have an analogous geometry as seen in Figures [1] and [2j The networks 
involving cycles of different types or of the same type but with a different 
number of nodes are represented in Figured 

4 Stability of cycles in type A networks 

In this section, we calculate the c-indices for each network of type A. We 
point out that, even though the geometry of the networks of type A and of 
type B is analogous, the stability indices behave differently. 

^We thank an anonymous referee for this comment on the superscripts for A cycles. 
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Figure 1: The cycles in the {A^jA^)- and the networks. The 

common connecting trajectory is in P 12 . 



Figure 2: The cycles in the {A^ , A^)- and the {B^ , B^ ) networks. The 
common connecting trajectory is in P 12 . 

For the rest of this paper, we assume that trajectories which leave a 
neighbourhood of the whole network do not come back. This ensures that 
the local and non-local stability indices coincide, as is implicit in [25]. 

In order to calculate the stability index, we linearize the vector held 
at each node. At a node ^i, there are four eigenvalues, often named —r* 
(radial), —Cij (contracting), 6*^ (expanding) and tu (transverse) for the role 
they play in the geometry of a cycle. We differ from this notation only for 
transverse eigenvalues tu that are expanding or contracting with respect to 
some other cycle in the networkH The constants r*, Cij and Cik are always 
assumed positive but tu can have either sign. For transverse eigenvalues tu 
with eigendirections away from the network, whenever possible we assume 
tu < ^ik so that, when the transverse eigenvalue is positive, it is weaker than 
the expanding one. Dehne, following [22], a, = Cjj/ejfc and 6* = —tu/eik- 

Since the calculations in this section rely heavily on Theorem 4.1 in [25] . 

®These are then called eu and —Cu, respectively. 
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Figure 3: The or (left) and the ,C^)- or 

(AgjAgjAj) networks (right). On the right, the cycle {C^) has two 
connecting trajectories in common with each of the A'^ cycles {B^). The 
connection [^i —>■ ^ 2 ] is common to all three cycles. 


we provide its statement together with relevant information for its under¬ 
standing next. 

We begin with the functions hij for 1 < j < m and I < j from [251 P- 
900]: 

hjjiy) = y, 

{ +CX) if ai-bi < 0, 

-.'.if.M-.+i if o<o,-6,<l,, 

aihi+ij{y) - bi if ai - k > 1. 

Note that as usual the indices are to be understood modulo m, which is the 
number of equilibria in the cycle. With this we can reproduce the result on 
the stability indices of type A cycles. 

Theorem (4.1 in [25]) For a type A cycle with m equilibria the stability 
indices are as follows: 

(a) // p > 1 and bj > 0 for all j, then aj = -|-cxd for all j. 

(b) If p > 1, bj > —1 for all j and bj < 0 for j = Ji,... Jl, then 

cr, = min h-, „ f —— | — 1, where ?' 

^ S=.h,...JL \ bsj 


J if j < s, 

j — m if j > s. 
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(c) If p < 1 or there exists j with bj < —1, then aj = —oo for all j. 

When the cycles are part of a network the next result shows how their 
indices are further constrained. The next four results are proved for a generic 
connection [,^j —)■ and their four relevant eigenvalues have j as a hrst 

subscript. We preserve the notation of |2S] for the stability indices. 

Theorem 4.1. Consider a network of type A cycles with linearization at 
as above. Let —)■ ,^j] be a common connecting trajectory. At least one of 

the cycles has all c-indices equal to —oo. 

Proof. Assume, without loss of generality, that at fj the expanding eigen¬ 
values are Cja and Cjb satisfying Cja > Cjb- Let Xa and Xb denote the cycles 
corresponding to the eigenvalues Cja and Cjb, respectively. 

It suffices to show that for cycle Xb, we have case (c) of Theorem 4.1 in 
[25] . hence aj = —oo. In fact, for the cycle Xb at fj, we have tjb = Cja so 
that bj = —Cja/Cjb < —1. □ 

Note that for networks consisting of three cycles, two of these cycles have 
c-indices all equal to — oo. 

On the side of stability, we have the following: 

Lemma 4.2. For an A cycle, all finite stability indices are non-negative. 

Proof. In Theorem 4.1. of [25] the only case with finite indices is (b), where 
aj - = 0 for all j and thus aj = aj^+ > 0 . □ 

Taking into account the fact that for any cycle in aj = —oo for some 
j if and only if aj = —oo for all j (Corollary 4.1 of [2S]), this gives 

Corollary 4.3. Generically, an A cycle with aj > —oo for some j, is e.a.s. 

A simple observation of the computations involved in the calculation of 
the images of the maps hij provides the following result. Having established 
that if stability indices are not —oo they will be positive, we next address 
the question of whether they are hnite or inhnite. 

Lemma 4.4. For a type A cycle in a simple network in 

(i) > 0 ^ aj < -|-oo, 

(a) if tk G (0, Cfc) V/c 7 ^ j, then: aj = -|-oo tj < —Cj. 
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Proof. (i) Let tj > 0. Then bj < 0, which excludes case (a) of Theorem 
4.1 in [23]. In case (c) all indices are equal to —cxo, so we look at case 
(b), where 



(ii) If tj < 0 is the only negative transverse eigenvalue, the inequality 
0 < tk < Ck for all k j puts us in case (b) again, where now 



For all to be equal to +oo, we need Oj — bj < 0, so Cj + tj <0. □ 

So understanding a network of A cycles means distinguishing e.a.s. cycles 
from those with all stability indices equal to —oo. Notice that if all con¬ 
nections in a cycle have stability indices equal to —oo then almost all initial 
conditions near this cycle will either approach another cycle in the network 
or leave any neighbourhood of the network. For an e.a.s. cycle, it is of in¬ 
terest to determine which indices are hnite and which are equal to -|-oo. To 
simplify notation, we implicitly assume that the stability index is hnite when 
writing cXj > 0. 

Following the terminology of [23] we use the notation 



Note that p is different for each of the cycles and p > 1 is necessary for the 
respective cycle to be anything other than completely unstable. 


4.1 The network 


Denote by —>■ the common connecting trajectory of the network. Since 
this connection is contained in a plane, the linearization at and may be 
written, respectively, as 



where all constants are positive. Denote the two cycles by and X 4 . With¬ 
out loss of generality, assume that at ^b the positive eigenvalues Cfes and cm 
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take trajectories to and X 4 , respectively. Then the transverse eigenvalnes 
for X 3 are — 0^4 and 6 ^ 4 , while those for X 4 are — 0^3 and e^s. The assnmption 
that et ,3 > eu implies that X 4 has all stability indices eqnal to —00. The 
network is the nnion of the cycles depicted in Fignre[TJ 
The necessary data for Theorem 4.1 in [25] is as follows: 

at Ca, CLa = — > 0 and ba = — > 0 , 
ea 2 ea 2 

at ^b, cib = — > 0 and bb = —— G (—1, 0). 

CbS 6^3 

Proposition 4.5. //p > 1, the c-indices for X 3 in the network 

are: 

I +00 ^ Ca3 < Ca4 

(^ab ^ U, (Tba 1 

^ 0 Ca3 > Cai 

Proof. There is only one negative transverse eigenvalne for X 3 , which is — 0^4 
at fa- The transverse eigenvalne at fb is and since we assnme Cbs > e^, 
we are in case (ii) of Lemma ITTl so aba = +00 is decided by Ca3 ^ 0^4. The 
remaining index is necessarily positive becanse it cannot be eqnal to — cxd. □ 


4.2 The network (^43 , ^43 ) 


The network is geometrically as that in Fignre [2J We then have a cycle 
containing the nodes .^ 1 , .^2 and .^3 with a common connection with another 
cycle with nodes .^ 1 , ^2 and .^ 4 . We distingnish between these cycles by 
referring to the node they do not have in common. With this convention we 
talk abont the .^3 cycle and the .^4 cycle. We assnme, as nsnal and withont 
loss of generality, that 623 > 624 . According to Theorem 14.11 at least one 
cycle has all c-indices eqnal to — 00 . Under the assnmption that 623 > 624 , 
this is the ^4 cycle. In what follows, we calcnlate the expressions for the 
c-indices of the ^3 cycle, nsing Theorem 4.1 in Podvigina and Ashwin Eg. 

The linearization at each node is the same as in [TS] (for a network of 
type B cycles). 


Proposition 4.6. Assume p > 1. The ^3 cycle is e.a.s. and its c-stability 
indices are as follows: 


(J12 > 0; (731 


-|-oo 


> 0 



Ci4 > Ci3 
Ci4 < Ci3 


<723 


-|-oo 


> 0 



(ci4 > Ci3 V C32 < C34) A C34 > 0 
(ci4 < Ci3 A C32 > C34) V C34 < 0 
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Proof. The required information at the nodes of the cycle is as follows. At 
.^1 we have 


At 6 , 


At 6 , 


Oi — 


and 

bi = 

^>0 


ei2 



ei2 

0.2 — 

C21 

and 

62 = 

624 ^ 


623 



623 

03 = 

C32 

and 

b3 = 

C34 


631 



631' 


The proof has to proceed in two cases, according to the sign of 63 . 


Case 63 > 0 : This occurs when C 34 > 0. We then apply Theorem 4.1 (b) 
in |25]. Along the connection [^3 —)• ^ 1 ] we have 

<^31 = ^ 1,2 

The value of hi ^2 depends on the sign of ai — 61 = If C 13 — C 14 < 0 

then hi ^2 = +00 and so is CT 31 . Otherwise, hi ^2 is hnite and a^i > 0. 

For the c-index along [^1 ^ 2 ] we have 



Finally, 

^■23 = ^3-3,2 ~ ^ ~ 

This is equal to +cxd if < 0, which is the same as C32 < C34. 

Otherwise ho ,2 is a function of hi, 2 , giving the same conditions as for < 731 . 




Case 63 < 0 : This occurs when C 34 < 0. We then apply Theorem 4.1 (b) 
in but now have to account for two negative 6 /s. In this case, is the 
minimum of two images of hij. Using Lemma 031 we know that all c-indices 
are positive and hnite, with the possible exception of cjjj where the transverse 
eigenvalue at is negative. This is possible only for CT 31 , so Lemma 031 yields 

0-31 = +CX) ci3 < C14. □ 

This is very different from the geometrically identical network 

for which various stability conhgurations (e.a.s. and non-e.a.s.) may occur, 
see section 5 in [ 6 ]. 
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6 cycle 

6 cycle 

A ^ cycle 


Ci3 ~ Ci4 

ai = —; bi = — 
ei2 ei2 

Cl4 , Ci3 

Oi = —; bi = — 

612 612 

Ci4 - 643 

oi = —; bi = — 
642 642 


C21 ~ 624 

02 — —; 02 —- 

623 623 

C21 , 623 

02 — —; 02 —- 

624 624 

C24 Y 624 

02 — —; 02 —- 

623 623 

^3 

C32 r 634 

03 = —; 63 =- 

631 631 


C32 r 634 

03 = —; 63 =- 

634 634 

6 


C42 , C43 

04 — —; 04 — — 
641 641 

C43 - 642 

04 — —; 04 — — 
641 644 


Table 1: The quantities aj and bj for cycles in the {A^ ,^ 43 ,^ 4 ) network. 


4.3 The network (^43 

Consider a network consisting of an cycle [.^1 —)■ .^2 —t ^3 ^4 —)■ .^ 1 ] and 
two A^ cycles, [.^1 ^2 ^ ^3 ^ ■Ci] and [.^1 —)■ ^2 ^ ■^4 .^ 1 ], which we 

call .^3 and .^4 cycle, respectively. Note that the geometry is not symmetric 
with respect to the A^ cycles. This is the only way to build such a network, 
because there are only four eigendirections at each equilibrium. 

We determine when each of the cycles is e.a.s. - for ease of reference we 
have listed the quantities aj and bj in Table [H We use aj, bj, dij ,... for the 
^3 cycle, plain letters for the ^4 cycle and aj, bj, aij, ... for the cycle. 

Lemma 4.7. In the {A'^,A^,A^) network 


(a) the ^3 cycle is e.a.s. if and only if p > 1 and 624 < 623 and 634 < 634 . 
If this is the case, then its stability indices are 


csi — 


+ CX) Ci 4 > Ci 3 

>0 Ci4 < C13 


0 ' 12 , <^23 > 0 . 

(b) the ^4 cycle is e.a.s. if and only if p > 1 and 624 > 623 . If this is the 
case, then its stability indices are 


0'4l5 <^24 — 


+ CX) Ci 3 > Ci 4 

> 0 C43 < C44 


cri2 > 0. 
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(c) the cycle is e.a.s. if and only if p > 1 and 624 < 623 and 634 > 631. 
If this is the case, then its stability indices are 


_ J 

f +CO 


Cl3 > Ci4 




[>0 


Cl3 < Ci4 



_ J 

f +00 


Cl3 > Ci4 

V 

642 > C43 

0-34 = < 

[>0 


Cl3 < Ci4 

A 

642 < C43 

<^12, <723 > 0. 






Proof. In all cases, the conditions for e.a.s. are straightforward applications of 
Theorem 4.1 in [ 25 ]: they are necessary and snfficient to have bj,bj,bj > — 1 
for all j in the respective cycle. So we calculate the indices for each case, 
making use of Table [TJ 

In case (a) we look at the .^3 cycle. Its only negative transverse eigenvalue 
is —Ci4 at .^1. So Lemma ITTI immediately yields a'i 2 , 5'23 > 0 , and dsi = +00 
if and only if C14 > C13. 

In case (b) the ^4 cycle has two negative transverse eigenvalues, thus 
62 < 0 is the only negative bj. From Theorem 4.1 in [ 25 ] we get 

1 , 

1 = £21 _ 1 > 0, 

623 

1 . 

Now hi 2(.) is equal to +cxo if ai — bi = < 0 . Otherwise it is a hnite 

expression involving h2p{.), which is hnite. For ho,2(-) we have 02 — 62 = 

^ ^ hnite expression involving hi^2(-) and thus the reasoning 

from above applies. 

For statement (c) and the Cf cycle we have 62, 63 < 0 and thus 
di 2 < hi ,2 (I - 1 = — - 1 < + 00 , 

d23 < h3,3 ( -i I - 1 = — - 1 < + 00 , 

V 63 / 631 

0-34 = min (ho,2(-l/62), ho,3(-l/63)) - 1, 
d4i = min (hi,2(-l/62), hi,3(-l/63)) - 1. 
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Conditions for the latter two indices to be finite follow in the same way as 
before. □ 

No two cycles can be e.a.s. simnltaneonsly. This is clear even withont 
lemma ITTI because they all share the trajectory [,^i —)■ ,^2]- We add to this 

Corollary 4.8. As long as p > 1 for all cycles, the (Ag , Af, Af) network is 
not completely unstable. 

The network Af): This network may be thought of as obtained from 
the previous one by deletion of one of the Af cycles. If we keep the assump¬ 
tions on the signs and magnitude of the transverse eigenvalues, the c-indices 
are the same as above. In this case, however, it is possible to make the net¬ 
work potentially more stable by choosing transverse eigenvalues to be neg¬ 
ative whenever possible. Assume we delete the ^4 cycle, i.e. the connection 
[^2 ^4], to obtain the network (Ag , Af). We may now choose ^2 = ^24 < 0 . 

Both remaining cycles are affected in the same way. For the ^3 cycle we 
may now have 02 — 62 < 0 (if and only if C21 < —624), then di2 = +00. 
Similarly, for the Af cycle we get di2 = +cxd if and only if C21 > —624. 

5 Concluding remarks 

We have contributed to the study of a certain type of simple heteroclinic 
networks in not just by providing a complete list of such simple net¬ 
works, but by studying the stability of cycles in networks of type A. In the 
heteroclinic context, type A networks have been overlooked. We show that 
even though their geometry resembles that of networks of other types, the 
stability of A cycles can be very different. 

The present list of simple heteroclinic networks provides convenient ex¬ 
amples for extending research on heteroclinic dynamics to questions arising 
from bifurcation, noise addition and so on. This is beyond the scope of the 
present paper. 
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A.l The network 

From the proof of Theorem 13.81 we know that P = (ki 2 ,/^i 3 ) — ^2 supports 
A 2 cycles. Since K 14 = ku o K 13 , there are three planar fixed-point subspaces 

P12 = {(a:i,a:2,0,0)} = Fix(Ki2) 

-P13 = {(xi,0,a;3,0)} = Fix(Ki3) 

.P 14 = {(xi,0,0,X4)} = Fix(Ki4), 

where (kia;) — Z 2 . The only fixed-point axis is 


Li = {(xi, 0, 0, 0)} = Fix(fi:i 2 , K 13 ) = Fix(P). 

We can produce a heteroclinic cycle of type in a similar way to that of 
the construction of the cycle in [ 6 ]. A P-equivariant vector field is 

4 

xi = aiXi ^ biix'^ + cixl 

i=l 


X2 




X 4 


Q‘2'^2 + 


03 X 3 -|- 


CI4X4 


hixj^ 


X2 C2X2X3X4 


X3 + C3X2X3X4 


X4 -F C4X2X3X4. 
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The coefficients ha need to be different from the rest bnt the remaining three 
may be eqnal|§ 

Choose tti > 0 so that the origin is nnstable. We can then choose bh — 
4aiCi > 0 as well as bu < 0 and ci < 0 to ensnre the existence of only two 
distinct eqnilibria on the xi-axis. The remaining coefficients may be chosen 
to provide the necessary saddle-sink connections. 

Therefore, this normal form snpports an ^4^ cycle and two of these can 
be put together in an network. 


A.2 The (AgjAg) and networks 

For the construction of a cycle of type we look at 

r = {ni 2 , ^ 13 ) 1^34), 


which contains all Kij and —Id G T, so that the same symmetry supports an 
cycle. A T-equivariant vector field for this symmetry is 


Xj — 


'^^bjixf I Xj + CjXiX2X^X4Xj. 


, i=l 


When restricted to a plane, say P 12 , we obtain 


xi = aiXi + I 'y'biixl 

K i=l 
^ 2 

X2 = a2X2 + I ^ b2ix‘^i 


Xi 


X2. 


s. i=l 


Note that this does not support a heteroclinic cycle of type as the two 
equilibria on the xi-axis would be related by symmetry: 

aiXi + biixl = 0 -v^ Xi = 0 or a:i = ±y/—Oi/fen, 

where r.(^-ai/6ii, 0, 0, 0) = {(0, 0, 0), 0, 0, 0)} is 

the corresponding group orbit. 

Coefficients can be chosen to ensure a saddle-sink connection in P 12 in a 
standard way. Two cycles of type Ag can be put together in a network in 
a way analogous to that used for the Pj cycle. Again the information at 
linear level required for the dynamics is the same as that obtained for the 
network. 

®This vector field is not equivariant for T = {k 2 , K 3 , K 4 ) used for the cycle (see [5], 
Appendix A). However, when restricted to the planes Pik the two vector fields do coincide. 
The information at linear level for the dynamics coincides with that for the cycle. 
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